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Within the QED effective action approach, we study 
the propagation of low-frequency hght at finite temperature. 
Starting from a general effective Lagrangian for slowly varying 
fields whose structure is solely dictated by Lorentz covariance 
and gauge invariance, we derive the light cone condition for 
light propagating in a thermalized QED vacuum. As an appli- 
cation, we calculate the velocity shifts, i.e., refractive indices 
of the vacuum, induced by thermalized fermions to one loop. 
We investigate various temperature domains and also include 
a background magnetic field. While low-temperature effects 
to one loop are exponentially damped by the electron mass, 
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there exists a maximum velocity shift of —Sv 
intermediate-temperature domain T ^ m. 
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I. INTRODUCTION 

As an application of the QED effective action approach 
at finite temperature (Q-Q, we investigate the aberration 
of low-frequency photons from the usual vacuum light 
cone in the presence of a heat bath. The propagation of 
photons can deviate from the classically free vacuum be- 
havior when external perturbations exert their influence 
on the vacuum polarization of the photon. In a per- 
turbative language, virtual processes, e.g., the electron- 
positron loop, confer the properties of the participant 
particles to the photon. By this means, the photon can 
therefore acquire, for instance, a "size" (of the order of 
the Compton wavelength of the participant particle) and 
a "charge distribution" . Finally, external perturbations, 
e.g., electromagnetic or gravitational fields, cavity- like 
boundary conditions, or a heat bath, can interact with 
these induced photon properties. The resulting modified 
photon propagation can be expressed in terms of a light 
cone condition which describes the deformation of the 
usual vacuum light cone as observed by the photons with 
frequency uj lower than the electron mass m. A variety 
of these effects have been under consideration in recent 
years (see, e.g., g and references cited therein). 

Concentrating on the deviation of the speed of light 
from the unmodified- vacuum velocity c (=1 in our 



units) by thermal effects, the velocity shiftN in the low- 
temperature region, T <C to, where m denotes the elec- 
tron mass, yields (5w^ ~ ~T^ /rn^. This important result 
can be obtained from the two-loop polarization tensor, 
in which the radiative photon within the fermion loop 
is considered to be thermalized 0i Hi- The same result 
has been found by an effective one-loop calculation em- 
ploying the Heisenberg-Euler Lagrangian p]n. The latter 
approach was also put forward in ||6|. Therein, the equa- 
tions of motion of the effective Lagrangian have been 
used to derive the light cone condition, and the effects 
of finite temperature have been introduced by evaluating 
thermal vacuum expectation values of the field quanti- 
ties. The high-temperature domain has also been stud- 
ied in Q by taking the effects of additionally thermalized 
fermions approximately into account. 

Since the effective action provides us with a convenient 
interface between the full quantum theory and classical 
field theory, it is the appropriate tool with which to de- 
scribe electromagnetic waves of low frequency {lo <C ttt-), 
which are the low-energy degrees of freedom of QED. 
The modified, e.g., thermalized, vacuum is viewed as 
a medium with optical properties which exert an influ- 
ence on the propagation of a plane wave. In the present 
work, we rigorously extend the effective action approach 
to the light cone condition as applied in j|] to the finite- 
temperature case. For this, the dependence of the ef- 
fective action on the complete set of gauge and Lorentz 
invariants of the given situation has to be considered; 
contrary to the zero-temperature case involving the two 
standard invariants T and Q, we have to take two fur- 
ther invariants into account which arise from an addi- 
tional four-vector associated with the heat bath. This is 
compulsory for an exact treatment of the problem in the 
effective action approach, and has not been considered in 
earlier works. The only input for deducing the light cone 
condition for low-frequency waves will be a general La- 
grangian whose dependence on electromagnetic field and 
temperature is solely dictated by Lorentz covariance and 
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^The velocity shift equals minus the shift of the refractive 
index associated with the modified vacuum. 

■^In this reference, the correct numerical result was found 
for the first time; furthermore, the similarities between the 
present problem and the Scharnhorst effect |lO|, i.e., light 
propagation in a Casimir vacuum, have been pointed out. 



gauge invariance. Therefore, the derivation of the hght 
cone condition does not rely on perturbation theory or a 
perturbative expansion of the Lagrangian. 

As an apphcation for our forniahsm, we fall back on 
perturbation theory and employ the thermal one-loop 
effective action of QED which allows for an investiga- 
tion of light propagation in a low- and an intermediate- 
temperature domain, T r^ m. It is a priori clear that the 
low-temperature velocity shift ^ T* will not be found to 
this order of calculation, since it is a two-loop effect in 
an exact field theoretic frameworkP . Unfortunately, a full 
two-loop calculation of the therma QED effective action 
has not yet been performed. At intermediate tempera- 
ture, T ^ m, the one-loop contribution is expected to be 
the dominant one. Apart from soft logarithmic correc- 
tions, we find a maximum velocity shift which reduces 
the hght velocity by 6vl^^^ = -^■ 

Our results do not apply to values of temperature 
at which the physics of the e"'"e~-gas is dominated by 
plasma effects, T ^ m. Then, signal propagation will 
completely be determined by plasma modes which are 
not considered in our approach. Especially, waves of low- 
frequency cease to propagate, since the appearance of a 
plasma mass serves as a cutoff for low frequencies. 



II. EFFECTIVE ACTION AT FINITE 
TEMPERATURE 

The low-energy effective action of QED for slowly vary- 
ing fields at finite temperature and zero density can only 
depend on a small set of invariants. In particular, a de- 
pendence of the effective action on the (infinitely many) 
invariants involving derivatives of the fields can be ne- 
glected if the typical scale associated with the variation 
of the fields in spacetime is much larger than the Comp- 
ton wavelength of the electron (Ac = 1/m). Concerning 
light propagation, this is a reasonable assumption as long 
as the frequency of the propagating wave is much smaller 
than the electron mass [lo <^ m) and possible additional 
background fields are almost homogeneous. 

In the sense of relativistic equilibrium thermodynam- 
ics, the Lorentz covariant and gauge invariant building 
blocks of these invariants are given by the field strength 
tensor and its dual, F'^", *F^^'' , and by the heat-bath vec- 
tor n^. The latter is on the one hand characterized by 
the value of its invariant scalar productFI, n^n^t = — T^, 
where T denotes the heat-bath temperature, and on the 



^In p], the loop counting refers to the order of the 
Heisenberg-Euler Lagrangian employed during the approxi- 
mation process for thermalizing the fermions and therefore 
differs from the present one. 

*We employ the metric g=diag(-l, 1,1,1). 



other hand related to the heat-bath 4- velocity u^^ via the 
invariant parameter T: n^ =Tu^. 

From these objects, we can construct the following set 
of invariants, which represents a conventional choice: 

T= V''F^.-i(B2-E2), 

g ^ l*F^'-F^, = -E • B, 

£ := ;^(naF"'')(n/3F^) ee {u^F<^^){upF^ ^) , 
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(1) 



Incidentally, parity invariance requires the action to be 
an even function of Q and prohibits the generation of a 
"Chern-Simons"-like term ~ u^*F^^'^ A^. Without loss of 
generality, we assume Q > Q for reasons of uniqueness. 
The convenience of the choice of £ becomes clear in the 
heat-bath rest frame, where we find f = E^. It is easy to 
check that the maximum number of classical invariants 
obeying these symmetry requirements is four (see, e.g. 
Sec. Ill of Ref. ||^); the set of Eq. (|l|) is thus complete. 
From these symmetry considerations, we conclude that 
the low-energy effective Lagrangian of QED for slowly 
varying fields at finite temperature can generally be writ- 
ten ag|^ 



C = C{T,g,£,T) 



(2) 



In other words, we are dealing with the local part (no 
derivatives of the fields) of the complete non-local effec- 
tive action (involving infinitely many derivative terms of 
the fields). Of course, even the local part of the exact ef- 
fective action incorporating the contributions of all loops 
is presently out of reach. Hence in practical applica- 
tions, we will rely on a perturbative expansion of Eq. 
(g); in fact, only the lowest order has been identified up 
to now. The derivation of this one-loop QED effective 
action for arbitrary constant electromagnetic fields at fi- 
nite temperature has required more effort than the purely 
magnetic case |^1|. The first comprehensive study has 
been elaborated by Elmfors and Skagerstam pj employ- 
ing the real-time formalism. For our purposes, we make 



^Apart from the implicit (and soft) spacetime dependence of 
the fields, the system under consideration is homogeneous in 
space and time (thermal equilibrium). Therefore, the effec- 
tive Lagrangian cannot depend explicitly on the coordinates 
x'^ . Although the finite-temperature case closely resembles 
the case of a Casimir vacuum, the latter is not invariant un- 
der translations orthogonal to the Casimir plates; hence, the 
effective Lagrangian for the Casimir vacuum can and, in fact, 
does depend explicitly on the spacetime coordinates llU]. In 
order to perform the considerations of the present work for 
the Ccisimir vacuum, we thus had to include another invari- 
ant of the form x'^riu where n'^ denotes the normal vector 
associated with the Casimir plates. 



use of the representation of the effective action as given 
in [g| . Therein, the field and temperature dependence is 
expressed in terms of the complete set of invariants as 
given in Eq. (|l|). 

Furthermore, introducing the abbreviations 



V-^' + 9^ + 3" 



b:= (/f^T^-^ 



(3) 



the thermal contribution to the one-loop effective La- 
grangian at zero density reads ||] 
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whereby the exponent is given by 
h[s) = .^ ,^ eacoteas 
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(5) 



The complete effective Lagrangian then consists of £ = 
Cci + ^^+^^'^+^j, where Cd = —J- denotes the classical 
Lagrangian and C^ represents the well-known one-loop 
part |13| at zero temperature, i.e., the Heisenberg-Euler- 
Schwinger Lagrangian. For later use, we also incorpo- 
rated the free photonic contribution C^ which is equal to 
minus the free energy density of a free photon gas (Stefan- 
Boltzmann law), C-,{T) — j^T*. Incidentally, only the 
zero-temperature one-loop part C^ is affected by renor- 
malization; for fixing the counter-terms, we require the 
complete Lagrangian to meet with the pure classical part 
in the weak- field limit (at T = 0). 



III. LIGHT CONE CONDITION 

The system under consideration is a propagating plane 
wave field in the presence of a vacuum that is modified by 
an external constant electromagnetic field at finite tem- 
perature and zero density. Following M, we assume the 
propagating wave to be of low frequency and neglect any 
vacuum modifications caused by the propagating light 
itself. In other words, we neglect dispersive effects and 
exclude non-linear self-interactions of the wave. Thus, 
we can linearize the field equations with respect to this 
plane wave field. 

Within the framework of these assumptions, the field 
equations are obtained from the Euler-Lagrange equa- 
tions of motion of the effective thermal QED Lagrangian: 
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(6) 



The Lagrangian C that we are going to insert into Eq. (g) 
contains the classical part as well as the zero-temperature 
and thermal one-loop contributions: £ = Cc\ + C^ + C^^; 
however, the following derivation of the light cone condi- 
tion does not rely on any perturbative approximation of 
C. It is only necessary that the field dependence of the 
Lagrangian be completely contained in a set of three lin- 
early independent invariants for which we take the stan- 
dard invariants J^, Q and the invariant £ as defined in Eq. 
(Q). Note that the differentiation with respect to F^i, has 
to be performed with regard to its anti-symmetry prop- 
erties. 

After moving the space-time derivative to the right em- 
ploying the Bianchi-identity 9^*F'"' = 0, we may write 
Eq. (^) in the form 



Q^d^Cd^F^'' 



d,F-P ( Im'^I 



(7) 



where M^'i denotes the tensor. 
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Faj} + 2dg£Cu[^(uF)p^ 



-Au^^'{uf) 



—^Fap H — Faf} +2d£Cuia{uf)i3] 



(8) 



where we employed matrix notation, (Ffc)"^ = F'^ak", or 
(ku) = k^Ua, and a prescription for anti-symmetrizing 
tensor indices: A^^"'^ = \{Ai"' - A"''). Note that the 
tensor M^^ is anti-symmetric with respect to its upper 
as well as its lower indices, and contains first and second 
derivatives of the Lagrangian with respect to the three 
invariants. As mentioned above, there is no explicit de- 
pendence of C on the coordinates because of homogeneity 
in spacetime. Equation (Q) represents the effective field 
equations for slowly varying fields in a thermalized QED 
vacuum. It replaces the linear Maxwell equations in vac- 
uum dfj,F'^'^ = 0, defining a new theory of "classical" 
non-linear electrodynamics. 

Following the techniques of Q, we adapt these equa- 
tions to the case of a propagating plane wave: first, we de- 
compose the electromagnetic field strength into the field 
strength of the plane wave field /'''' and an additional 
homogeneous background field F^'^: Ff"" -^ F^'' + f" . 
Now derivatives of the field act only on the plane wave 
part: (9^(F"^ + f"'^) = df,f"^. Transition to Fourier 
space leads us to the replacement —id^f"^ -^ k^f"^^ ^ 
where k^ denotes the wave vector of the plane wave field 
/"^. The latter is proportional to f"!^ ~ fc"e^ - fc^e", 
whereby e" represents a polarization vector of the plane 
wave. Furthermore, we adopt the Lorentz gauge condi- 
tion for the plane wave field and sum over polarization 



states; the latter step is, on the one hand, convenient, 
and, on the other hand, reasonable for a purely thermal- 
ized vacuum which is naturally isotropic. In the case 
of an additional background field, we lose information 
about effects of birefringence by summing over polariza- 
tion states (cf. later). Finally, we arrive at the desired 
light cone condition: 



= 2 d:FC fc2 + 2 
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(9) 



Since we have linearized with respect to the propagat- 
ing plane wave field, the field quantities appearing in Eq. 
(^ now describe an externally applied background field 
only. Following summation over polarization states, the 
field strength tensor of the plane wave field has com- 
pletely dropped out of the formula. This corresponds to 
our approximation of neglecting vacuum modifications 
caused by the plane wave field itself. In this sense, the 
plane wave field resembles a small test charge of classical 
electrodynamics. 

Evaluating the contractions of Af^^ with fc^/c", we find 
the light cone condition of a soft plane wave propagating 
in an external field at finite temperature in its final form, 

0- (a^/:+ gd:FgC-Tdlc)k'+ i(a|.+ dl)c{FkY{?k), 

-d£Ck'^ + 2[d£C + £dlC]{kuf (10) 

+2dr£C {ku){VkY {?u)u + 2dg£C (fcu)(*FA:)"(FM), 
+2{d:F£C - dlC) {uVkf + 2dg£C (uF/c)(u*FA:). 

Let us stress once more that the field invariants !F,g,£ 
in this equation characterize an externally applied back- 
ground field only. The first line in Eq. (Iln) corresponds 
to the purely field-modified light cone condition as de- 
rived in B , which can conveniently be expressed in terms 
of the energy-momentum tensor (or its vacuum expecta- 
tion value) of the electromagnetic field. 

In the following, we simply confine ourselves to the 
case of vanishing external field: F'^'^,J^,Q,£ — > 0. In 
this limit, the light cone condition becomes extremely 
simplified, and, except for the term linear in J- which is 
filtered out by i9j^>C|jp-,£;.£=o, only the first term in square 
brackets in Eq. ( p^ survives: 

= {djrC-d£C)k^ + 2d£C{kuf, for T,g,£~^0. 

(11) 

Introducing the phase velocity v ^ t^^, which is identi- 
cal to the group velocity in the low- frequency limit, we 
may rewrite the light cone condition in the heat-bath rest 
frame ((fcu)^ — (fc°)^) in terms of the squared velocity, 



v^ = 



1 I 2d£C 



(12) 



In the present case of light propagation in a thermal- 
ized vacuum, it is physically reasonable to expect that 



V < 1 holds rigorously, since the propagating photons do 
not only virtually interact with the usual vacuum modes 
(zero-point fluctuations) but also with the thermally ex- 
cited modes. This may be interpreted as a kind of "resis- 
tance" implying a decrease of the propagation velocity. 
In order to maintain w < 1, the fraction in the denomina- 
tor of Eq. ( p^ ) should always be positive. Note that this 
statement is independent of any loop approximation. 

Now we turn to perturbation theory, where C can be 
assumed to consist of the classical term Cd = —T plus 
quantum correction terms such as those discussed in Sec. 
II. The latter are considered to be small compared to 
the classical one; hence, Eq. (p2|) to lowest order in the 
perturbation simplifies to: 



v^ ■^\-2d£li. 



(13) 



As will be demonstrated below, this formula describes 
sufficiently the thermally induced velocity shift to one 
loop for reasonable values of temperature. 

For the remainder of the section, we will transcribe 
the light cone condition Eq. (^l|) into a form which has 
proved useful for gaining intuitive insight into the prob- 
lem of light propagation in modified vacua [p|. For this, 
we introduce the vacuum expectation value of the energy- 
momentum tensor, 



l(r^^) =Lg^'' ^TdrCu^v!', 



(14) 



which can be obtained by varying the effective action 
with respect to the metric g^j/, and treating the vari- 
able r as T = ^—n^^gy^yU^ . Equation (nj) is the finite- 
temperature analogue of Eq. (19) of H where the vacuum 
expectation value of T'^'^ for an electromagnetic back- 
ground has been considered. Solving Eq. (nj) for u'^u'^ , 
we may substitute this into Eq. (hlf) and arrive at: 



(15) 
Qt for the 



fc2 == Qt (T^ fc^fc,, 

where the so-called effective action charge 
present system is given byQ 

2 dsC 



-djrC + d£C+2,^C 



(16) 



From Eq. ( [Tsl ) , we can read off that the deformation of 
the light cone is determined by the (renormalized) energy 
density of the modified vacuum. The proportionality fac- 
tor Qt depends on the parameters of the system via the 



In H, the name "effective action charge" derives from the 
resemblance between the corresponding equation to Eq. ( |lq ) 
for electromagnetic fields and the Poisson equation, 2Q :— 
{d% + dg)C. This accidental feature does not hold in the 
present case; however, we stick to this nomenclature, since 
essential properties as well as the underlying physical picture 
also hold here. 



effective Lagrangian. Under the reasonable assumption 
that the deformation of the hght cone is bounded even 
in the high-temperature domain, one can conclude that 
the factor Qt has to decrease sufficiently fast for increas- 
ing temperature. Therefore, we expect Qt to be locally 
(charge-like) distributed in the parameter space of tem- 
perature. This justifies the nomenclature effective action 
charge. In fact, the one-loop results will confirm this 
picture. 



IV. ONE-LOOP RESULTS 



where the numbers stem from pure integrals over ana- 
lytic functions. Incidentally, the expansion for T ^ m 
is formally identical to Eq. (EO) with the factor of 
0.814. . . replaced by 1 as found in |lj]. 

For our purpose of investigating low- and intermedi- 
ate temperature domains, we observe that dj^C^^ {T) — 
O^aj-K) <C 1. For the calculation of the velocity shifts to 
order a/vr, it is thus sufficient to employ Eq. (|l3|), i.e., 
neglecting dj^C^'^ compared to dj^Cd = — 1. 

Turning to the low-temperature expansion of Eq. jl^), 
we have to take into account only the first term of the 
sum: 



As an application, we insert the one-loop effective ac- 
tion into the light cone condition. Let us first consider 
a thermalized vacuum with vanishing background fields. 
Therefore, the velocity shift stems from the thermal con- 
tribution C^^ only, and we need dsC = dgC^^ in the 
zero-field limit: 



dsC 
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(17) 



Regarding Eq. (O), we also need dj^C; for this, note 
that the zero-temperature one-loop part is rcnormalizcd 
in such a way that the term linear in T vanishes in 
order to recover the pure Maxwell theory in the weak- 
field limit. Hence, it is only the thermal contribution 
C^^ which provides for an additional term linear in T: 
djrL = —1 -I- dyrC^^ J where the (—1) stems from the 
Maxwell part Cd- In the zero-field limit, we get: 



9^£i^ = -|^^(-l)"Xo(-„). 

n— 1 



(18) 



For Eqs. (17) and (|18|), we employed the fact that sum- 
mation and proper-time integration in Eq. (0) are in- 
terchangeable in the zero-field limit. The proper-time 
integrals can then be identified as a representation of the 
modified Bessel function Ki and Kq fl^ . 

In the limiting cases of low, intermediate and high tem- 
perature, Eqs. ( |l7|) and (|l8| ) can be expanded in terms 
of the parameter — . Let us first consider Eq. ([l8|): at 
low temperature, only the first term of the sum must be 
taken into account, which leads to 



drC^^{T^m)~ — J e- 

Stt \ m 



(19) 



For an intermediate- or high-temperature expansion, the 
infinite sum must be completely taken into account; this 
can be achieved employing the techniques proposed in 
the appendix B of |||. The result for Eq. (18) at T ~ m 



d,C'^(T«m)^i-J-J^e- 



a 
3^ 



(21) 



Inserting this into Eq. (13), wc find for the squared veloc- 
ity at low temperature an exponentially decreasing mod- 
ification: 



v'^{T ^m) ~ 1 



TT y T 



(22) 



The high-temperature expansion of Eq. (|l7| ) can be 
worked out along the lines of App. B of [p[, leading to 



dsC' (T > m) 



(l-fc2^+0(f^)), (23) 



where /c2 denotes the number k2 = 0.213139 .. . and 
arises from a parameter-independent integral over an- 
alytic functions during the expansion. The light cone 
condition then yields for the squared velocities at high 
temperature: 



v^{T-»m) ~ 1 



a 
3^ 



1 



' y2 



(24) 
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which implies a maximum velocity shift of Sv 

a ^ 1 

3tt — 1291 ■ 

It should, however, be noted that an expansion for 
T/m ^ 1 is a formal trick to extract analytical results. 
For identifying the values of temperature to which the 
light cone condition Eq. ( [lO| ) is applicable, we notice 
that the frequency of the plane wave should, on the one 
hand, be smaller than the electron mass in order to justify 
the assumption of slowly varying fields, and on the other 
hand, be larger than the plasma frequency in order to 
ensure the existence of such a propagating mode: 

Wp < w < TO. (25) 

For the plasma frequency corresponding to the Debye 
screening mass, we employ the representation found in 

i: 



d^C'^ = =^ 



2a 
3^ 



(0.666...) + (0.814... ) In — 

ra 



(20) 



w; 



-TO^5^(-l)"X2(fn). 



(26) 



The maximum value of temperature up to which the 
low-frequency assumption can formally be justified is 
determined by u;p(Tmax) = rn. Numerically, one finds 
Tma.x/'m- — 5.74 .... Of course, this is just a formal value; 
in order to obtain reasonable results, i.e., in order to sat- 
isfy relation (|2^) , the actual temperature should be kept 
smaller than this maximum value. 

Nevertheless, the numerical results given below con- 
firm that the formal expansion for T/m ^ 1 is already 
appropriate for T ^ m, which justifies extracting physi- 
cal conclusions from this analytical procedure. 

For the remainder of the section, we additionally al- 
low for a further perturbation of the vacuum in form 

2 

of an external weak magnetic field, B ^ B^- = 2i_, j^ 
terms of the invariants, purely magnetic fields imply that 
G,£ -^ 0,JF^ ^^^: where the latter expression is under- 
stood to be valid in the heat-bath rest frame. Due to the 
simple form of the four-velocity vector u^ — (1, 0, 0, 0) in 
this special frame, terms containing F^^u^, vanish, since 
this product is proportional to the electric field, which is 
assumed to be zero. 

Taking these considerations into account, the light 
cone condition (nfl) is reduced to: 

+2d£C{kuf. (27) 



These terms can be inserted into Eq. ( |27| ) together with 
the classical, the one-loop, and the purely thermal part 
of the effective Lagrangian, CchC^^ and C^'^ . Introduc- 
ing the phase velocity i" = -fej again and omitting the 
terms of higher order in T /m and eB /m?' , we arrive at 
the following polarization sum rule describing the mod- 
ification of light propagation in a QED vacuum at low 
temperature and weak magnetic background field: 




2tt m? 4 m^ 



,-f Q-2 



(30) 



In the first line of this equation, we encounter the famil- 
iar classical, purely thermal, and purely magnetic con- 
tribution. The second and third line describe the cor- 
rection terms due to the simultaneous presence of a heat 
bath and a magnetic background field, indicating that the 
combination of both perturbations renders more than a 
simple sum of the single effects. 



We will restrict the following investigations to experimen- 
tally more accessible parameter values, i.e., weak fields 
and low temperature (compared to the electron mass). 
Here we would like to answer the question of whether 
there is an additionally induced velocity shift for the 
case of combined magnetic and thermal vacuum mod- 
ifications. For each single case, the velocity shifts are 
known; e.g., the polarization-summed velocity shift of 
light propagating in a magnetic background field is given 
by [pi: 



45 TO* 



(28) 



where Q denotes the angle between the propagation di- 
rection and the magnetic field. The velocity shift for the 
purely thermal part is found in Eq. (2^). 

In order to find a contribution to the velocity shift for 
the combined case apart from the trivial sum of the single 
cases, we need an appropriate expansion of the thermal 
one-loop Lagrangian £^"^ for weak fields and low temper- 
ature. This expansion has been obtained in a systematic 
way in |I4|] ; the same result can be found by a direct 
Taylor expansion of Eq. (0). In the desired limit, the 
required terms which are quadratic in the field invariants 
and dominant at low temperatures read: 



ACbt{J'.Q,£,T) 



a 
45 



^T8 - 8J^2 _ 13^2 






T3 



(29) 



V. DISCUSSION 

A. Low^ Temperature 

At temperatures well below the fermion mass, the one- 
loop modification of the velocity of light as described by 
Eq. ( |2^ ) vanishes exponentially. It is the mass of the 
fermions in the loop that is responsible for the attenua- 
tion of thermal effects, since it counteracts thermal exci- 
tations of higher virtual modes. As a consequence, the 
usual hierarchy of loop calculations is inverted: the dom- 
inant thermal contribution to the low-temperature veloc- 
ity shift stems from the two-loop graph in which the inter- 
nal photon line is thermalized. Since the photon is mass- 
less, thermal excitations produce a low-temperature ve- 
locity shift of <5w = -^M^a^i^ §- i. Incidentally, the 



three-loop diagram involving two radiative thermalized 
photons of course contributes subdominantly ^ o?^-^, 
as can be deduced from the findings of [|l6| . 

However, we would like to point out that the dom- 
inance of the two-loop contribution docs not seem to 
hold over the complete range where T < tti, as can 
be discovered numerically. The increasing influence of 
the one-loop term is due to the factor of T^^l"^ and, 
of course, the proportionality to a (and not a^) in Eq. 
(p3). The one-loop contribution becomes dominant for 
comparably low temperatures of T jm > 0.058 (Fig. 1). 
Of course, this statement should be handled with care 
because here we are comparing a one-loop result with 



thermalized fermions to a two-loop results without ther- 
malized fermions. Nevertheless, an increasing two-loop 
contribution from thermalized fermions will always be 
suppressed by a factor of a, which could only be com- 
pensated for by an unexpected conspiration of numerical 
pre- factors. 
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FIG. 1. Low-temperature velocity shift 5v (=i5w^) in units 
of the vacuum velocity c = 1 versus the dimensionless tem- 
perature scale T /m\ the one-loop contribution as given in 
Eq. (E3) exceeds the well-known two-loop result ~ T^ jrrv' 
for comparably low values of temperature. 



B. High Temperature 

As discussed at the end of Sec. IV, the results of 
the formal high-temperature expansion should be ap- 
plied to values of temperature of the order of the elec- 
tron mass, T ^ m, where plasma effects ||lj do not yet 
dominate the physical properties of the modified vacuum. 
At even higher temperatures, T ^ m, the magnitude of 
the plasma frequency ujp serves as a cutoff for the low- 
frequency waves which belong to the main ingredients 
of our formalism. Then our calculations become mean- 
ingless, because the low-frequency modes simply do not 
propagate. 

Nevertheless, already in the temperature region where 
T ~ TO, our results indicate that the thermal excitations 
are no longer seriously damped by the mass of the fermion 
in the loop. The one-loop contribution obviously domi- 
nates the velocity shift; therefore, the usual loop hierar- 
chy is completely restored. 

It is remarkable that a maximally possible velocity 
shift to this order of calculation exists which is sim- 
ply given by i^w^^x = "^ (Fig- 2). Of course, this 
maximal velocity shift is only reached asymptotically, 
and therefore, strictly speaking, lies beyond the scope 
of the present formalism; nevertheless, the actual veloc- 
ity shift already comes close to the maximum value in 



the intermediate-temperature domain where T ~ to (Fig. 
2). To complete the discussion, it should be mentioned 
that the inclusion of the contributions from Eq. (|lq) in- 
crease the negative velocity shift proportional to % In — 
for T > 771. For reasonable values of temperature, this 
contribution is in fact completely negligible. 
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FIG. 2. Numerical evaluation of Eq. (113) in one-loop 
approximation: the thermally induced velocity shift —5v 
(= — ^Sv'^) in units of the vacuum velocity c = 1 is plot- 
ted versus the dimensionless temperature scale T/m; the left 
part of the curve is asymptotically described by Eq. (p2|), 
the right part by Eq. (M). The maximum velocity shift 
—5v = ^ (equivalently: —Sv = ^) is already approached 
at comparably low temperatures. 



C. Low Temperature and Weak Magnetic Field 

In order to justify the low-temperature expansion lead- 
ing to Eq. (p9|), we restrict the following discussion of Eq. 
( pO| ) to maximal values of temperature of Tmax/"^ — 1/3. 

Let us stress once more that the present results derive 
from a light cone condition which represents a sum-rule 
for different polarization modes. In fact, the QED vac- 
uum under the influence of a oriented background mag- 
netic field is reminiscent of an anisotropic medium with 
different refractive indices, i.e, velocity shifts, for each po- 
larization mode. In this sense, the velocity shifts found 
in Eq. ( pO| ) represent the arithmetic averages of those for 
the single polarization modes. A detailed discussion of 
magnetically induced birefringence at finite-temperature 
for which the eigenvalue problem of the tensor M^o (Eq. 
(P)) has to be solved is out of the scope of the present 
work. 

Concentrating on those terms in Eq. (pG) which de- 
pend on the magnetic field, we encounter the well-known 
zero-tcmpcrature velocity shift (last term of the first line) 



with its typical dependence on the angle Q =<)L(k, B); fur- 
thermore, there is also one thermal correction with this 
9-dependence (second line), and another which is inde- 
pendent of the direction of the magnetic field (third line). 
Although both thermal corrections vanish exponentially 
in the zero-temperature limit, the factors of T in the nu- 
merator provide for a strong increase for T/m > 0.1. 

Numerical analysis shows that the second line of Eq. 
( ^ ) contributes with an opposite sign compared to the 
zero-temperature part. Since both exhibit the same de- 
pendence on the magnetic field and the angle 0, we may 
add them and find that, on the one hand, thermal ef- 
fects diminish the zero-temperature velocity shift in a 
magnetic background by an amount of, e.g., ~20% at 
T/m ~0.25. On the other hand, the third line of Eq. (JS^) 
contributes with the same sign as the zero-temperature 
velocity shift, and becomes comparable to the latter (for 
Q = tt/2), e.g., ~ 80% of the zero-temperature velocity 
shift at T/m ~ 0.25. Therefore, at these values of tem- 
perature, the velocity shift loses its strong dependence on 
the direction of the magnetic field and the propagation 
direction and is partly dependent on the energy density 
of the magnetic field. 



D. Effective Action Charge 

In Eq. dl^), we were able to formulate the thermally 
induced deformation of the light cone in terms of the 
vacuum expectation value of the energy-momentum ten- 
sor times a factor Qt called effective action charge Q. 
The charge concept has been introduced as a useful pic- 
ture which provides for an intuitive understanding of this 
proportionality factor. Since any perturbed QED vac- 
uum can be expected to control the magnitude of the 
velocity shifts even for vacuum modifications of high en- 
ergy density, this factor Qt has to decrease sufBciently 
fast for increasing energy scale parameters (in this case: 
temperature). Therefore, the factor Qt should exhibit a 
localized distribution in this parameter space. 

For calculating the effective action charge Qt accord- 
ing to Eq. (y^ for a thermalized QED vacuum, we do not 
only consider the one-loop contribution from the thermal- 
ized fermions C^'^ , but also take the free photonic part 



An' 



TdTC-,{T) = --T 



L~f — fg-r* into account. Although it does not exert an 
influence on the velocity shift and drops out of Eq. (15), 



we have to include it in the present considerations in or- 
der to work with the complete vacuum expectation value 
of the energy- momentum tensor for the thermalized QED 
vacuum. 

To lowest order in a/7r, the formula for the effective 
action charge Eq. (p^) reduces to: 



Qi 



2dEL 



T OtC TOtC-, + TOtC^^ '' 



(31) 



where dgC is given in Eq. (17h, and the terms in the 
denominator can be written as: 



45 



TOtC' (T) 



T 



E 

n=l 



(-1)" 



Ksi^n 



(32) 



Note that the appearance of OtCj and OtC^^ in the de- 
nominator of Eq. ( |3l| ) corresponds to the appearance of 
a photonic and a fermionic part in the energy-momentum 
tensor. One usually expects the fermionic part to become 
important only for high temperature, T 3> m, where the 
fermions become ultra-relativistic. 
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FIG. 3. One-loop contribution to the effective action 
charge Qt in units of -^\ for high temperature, Qt decreases 
proportional to ^^ . The solid line corresponds to Qt as given 
in Eq. (pll); for the dashed line, the fermionic contributions 
^ TOtC have been omitted. In the low-temperature limit, 
the effective action charge to one loop vanishes due to the 
influence of a finite electron mass. 



On the one hand, we indeed find the expected localized 
behavior, as can be seen in Fig. 3: the effective action 
charge vanishes for high as well as low temperatures; in 
between, it develops a maximum at T/m ~ 0.22. On the 
other hand, it is interesting to note that the inclusion 
of the fermionic contributions ^ TOtC^'^ (solid line) in 
the denominator of Eq. ( p3l| ) becomes important for val- 
ues of temperature close to the maximum of Qt- This 
again indicates that the thermalization of the fermions 
becomes important already for comparably low values of 
temperature. 

Contrary to the cases discussed in [Q , the one- loop ef- 
fective action charge is not centered at the origin. This is 
because the electron mass damps the fermionic thermal 
fluctuations exponentially for small values of tempera- 
ture. 



VI. CONCLUDING REMARKS 

First, we would like to stress that the one- loop contri- 
butions to the velocity shift as calculated in the present 
work do not fit into the scheme of the "unified formula" 
proposed in ||] ; the latter connects the velocity shift with 
a shift of the vacuum energy density caused by external 
influences with a universal constant coefficient as propor- 
tionality factor. In the language of Eq. (^5|), this coeffi- 
cient has to be identified with the effective action charge 
Qt which in the present case is not constant at all but 
almost carries the complete physical information of the 
problem. This misfit indicates that the "unified formula" 
is an artefact of an approximation scheme rather than a 
fundamental principle |^. 

Furthermore, it should be pointed out that the max- 
imum velocity shift, — (Jw^ax — ^' cannot be viewed as 
an experimentally significant limiting value, since light 
propagation will successively be dominated by plasma ef- 
fects for increasing temperature aXT ^ m. Nevertheless, 
the existence and amount of such a maximum shift are at 
least interesting from a theoretical viewpoint, since they 
characterize the classically forbidden interaction between 
the modified vacuum and a photon which is exposed to 
the effects of vacuum polarization. Supposing there were 
no collective excitations constituting a plasmon for high 
temperature, then QED would not allow for an arbitrar- 
ily strong influence of a heat bath on the propagation of 
light for reasonable values of temperature. Of course, for 
extremely high temperature, the logarithmic corrections 
from dj^C in Eq. (|l2|) would also significantly increase 
the negative velocity shift, slowing down the speed of 
light; e.g., another shift of a/Zii would be reached at 
T/m ~ 6.4 • 10^. 

The problem of calculating the intermediate- /high- 
temperature velocity shift was also approached in [pj em- 
ploying a method reminiscent of the Born-Oppenheimer 
approximation: the coupling of the external field to the 
fermion loop was considered as an adiabatic interaction, 
while the photonic fluctuations were treated as the "fast" 
degrees of freedom. This calculation can be rated as 
a two-loop calculation, since the thermalization of the 
fermions has been approximately taken into account by 
using a thermal one-loop Lagrangian; photonic fluctu- 
ations were introduced in the same way as in the low- 
temperature case, i.e., by taking thermal vacuum ex- 
pectation values of the field quantities. Similarly to the 
present work, a maximum velocity shift has been found 
which is about three orders of magnitude smaller than 
the one-loop result given above. Indeed, this fits into the 
usual hierarchy of perturbation theory. 

To overcome the approximation scheme of B, which 
treats the thermalization of fermions and photons on a 
different footing, the present paper offers an appropri- 
ate tool by working with the complete set of invariants 
of the thermalized QED vacuum. However, for a direct 
verification of the results of B] for high temperature, a 



first-principles two- loop calculation of the effective action 
should be carried out in the same sense as suggested in 
the present work. 

Beyond the clarification of these technical questions, 
a two-loop calculation of the thermal effective action 
will be useful for studying further phenomenological as- 
pects of the physics of strong fields, such as thermally 
induced pair production or photon splitting. In view 
of the present results, one might speculate that, for 
these effects, two-loop contributions will be dominant at 
low temperature, since thermal one-loop corrections have 
proved to be vanishing (for pair production pUqJq] ) or of 
minor importance (for photon splitting U%). 
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